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Abstract 



Form factor axioms are derived in two dimensional integrable defect theories for matrix 
elements of operators localized both in the bulk and on the defect. The form factors of bulk 
53^1 operators are expressed in terms of the bulk form factors and the transmission factor. The 

, ■ structure of the form factors of defect operators is established in general, and explicitly 

calculated in particular, for the free boson and for some operator of the Lee- Yang model. 
^ _ Fusion method is also presented to generate boundary form factor solutions for a fused 

Q ■ boundary from the known unfused ones. 

o ■ 

^ ■ 1 Introduction 
O 

Q\ , The bootstrap program aims to classify and explictly solve 1+1 dimensional integrable quantum 
; field theories by constructing all of their Wightman functions (see [3] for a recent revievi^ and 
references therein). In the first step, called the S-matrix bootstrap, the scattering matrix, 
^ . connecting asymptotic in and out states, is determined from its properties such as factorizability, 
unitarity, crossing symmetry and Yang-Baxter equation (YBE) supplemented by the maximal 
analyticity assumption p]. In the second step, called the form factor bootstrap, matrix elements 
of local operators betv^^een asymptotic states are computed using their analyticity properties 
originating from the already computed S-matrix. Supposing maximal analyticity leads to a set 
of solutions each of v^^hich corresponds to a local operator of the theory [T]. In the third step 
these bulk form factors are used to build up the correlation (Wightman) functions via their 
spectral representations and describe the theory completely off mass shell. This program has 
been implemented for a wide range of theories (see [3]). . 

The analogous bootstrap program for 1+1 dimensional integrable boundary quantum field 
theories has been already developed. The first step is called the R-matrix bootstrap [5]: In 
boundary theories the asymptotic states are connected by the refiection (R)-matrix, which, as 
a consequence of integrability, factorizes and satisfies the unitarity, boundary crossing unitarity 
and the boundary YBE (BYBE) requirements. These equations supplemented by the maximal 
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analyticity assumptions make possible to determine the reflection matrices and provide the 
complete information about the theory on mass shell. In the second step we are interested in 
the matrix elements of local operators localized both in the bulk and also at the boundary. 
Due to the absence of translational invariance the bulk operators^ one point functions acquire 
nontrivial space dependence which can be calculated in the crossed channel using the knowledge 
of the boundary state together with the bulk form factors [22]. For the matrix elements of local 
boundary operators axioms can be derived from their analytical properties originating from 
the already computed R-matrix [9]. Supposing maximal analyticity leads to a set of solutions 
each of which corresponds to a local boundary operator of the theory and is uniquely related 
to a vector in the ultraviolet Hilbert space. The explicit form of the boundary form factors 
determine the boundary correlation functions via their spectral representation. This provides a 
partial description of the theory off the mass shell as a full description would include correlation 
functions of operators localized in the bulk as well, but this complicated problem has not been 
addressed yet. 

Since any two dimensional defect theory can be mapped to a boundary theory [13] the de- 
velopment of a separate bootstrap program for their solution seems to be redundant. However, 
integrable defect theories are severely restricted and one can go much beyond the boundary 
bootstrap program explained above: We can determine the form factors of both types of op- 
erators, those localized in the bulk and also the ones localized on the defect. With the help 
of these form factors we are able to derive spectral representation for any correlation function 
and in principle fully solve the theory off the mass shell as we will show in the present paper. 

In developing a defect form factor progam the first step is the T-matrix bootstrap. In- 
teracting integrable defect theories are purely transmitting [TOl [IH [12] and topological. As 
a consequence a momentum like quantity is conserved [HI [I5| and the location of the defect 
can be changed without affecting the spectrum of the theory [THIlT]. This fact, together with 
integrability leads to the factorization of scattering amplitudes into the product of pairwise 
scatterings and individual transmissions and enables one to determine the transmission factors 
from defect YBE (DYBE), unitarity and defect crossing unitarity [HI EHl [IS] . The second step 
is the defect form factor bootstrap: Once the transmission factors are known we can formulate 
the axioms that have to be satisfied by the matrix elements of local defect operators. We will 
analyze operators localized both in the bulk and also on the defect. By finding their solutions 
the spectral representation of any correlator can be determined and the theory can be solved 
completely. For simplicity we restricted our interest in the paper for theories with a single 
particle type. 

The paper is organized as follows: In section 2 we introduce asymptotic states in defect 
theories and the notion of the transmission matrix. Then the coordinate dependence of defect 
form factors is determined. By specifying the boundary form factor axioms we postulate the 
axioms for diagonal defect theories. In section 3 we determine the form factors of any operator 
localized in the bulk in terms of the transmission factor and the bulk form factors. For operators 
localized on the defect a procedure to calculate the general form factor solution is outlined. In 
section 4 we apply this technology to determine the defect form factors of the free boson and 
the Lee- Yang model. Finally we present how the fusion method can be adapted to generate 
boundary form factor solutions for a fused boundary from the known unfused ones. 
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2 Defect form factors 



In this section we present the axioms for the matrix elements of local operators between asymp- 
totic states. To shorthen the discussion we introduce Zamolodchikov-Faddeev (ZF) operators 
in order to describe both the multiparticle transmission process as well as the properties of the 
form factors. 

2.1 Asymptotic states and transmission matrix 

Multi-particle asymptotic states in integrable bulk theories can be formulated in terms of the 
ZF operators as 

|^i,...,^„) = A+(^i)...A+(e„)|0) 

All particles have different momenta Pi = msinh^j, thus in the remote past they are not 
interacting and form an initial state 6i > ... > 6'„. When time evolves they approach each 
other and after the consequtive scatterings they rearrange themselves into the opposite order: 

l^'l, ■ ■ ■ ,0n) =Yl^i^i~ ^j)\^n, ■ ■ ■ 

Here S{9) is the two particle scattering matrix which satisfies unitarity and crossing symmetry 

s{-e) = S'\e) ■ s{i7T -e) = s{e) 

This multi-particle scattering process easily formulated with the ZF algebra: 

A+{ei)A+{e2) = s{e, - ^2)a+(^2)^+(^i) + 27r6{ei -e^- tn) (i) 

where we extended their definition for imaginary 6 by postulating the crossing property 

A{e) = A+ie + tTc) (2) 

see [9] for the details. 

Once defects are introduced we have to make a distinction whether the particle arrives from 
the left (A) or from the right (B) to the defect. These particles can be even different from each 
other as they live in different subsystems. A multiparticle state is then described by 

1^1, ... , On, On+i, ...,e^) = A-^{e^) . . . . . . fi+(ft„)|o) 

where the ZF operators create particles on the right of the defect and satisfy similar defining 
relations to ^ with a possibly different scattering matrix. Yet, for simplicity, we restrict our 
discussion to the case when the two subsystems are identical with the same scattering matrix. 
Observe however, that this does not imply space parity invariance, since the defect may break 
it. In the initial state rapidities are ordered as 6*1 > . . . > ^„ > > 6n+i > . . . > Om- The final 
state, in which all scatterings and transmissions are already terminated, can be expressed in 
terms of the initial state via the multiparticle transmission matrix. 

n m 

|6'i, . . . ,9n] On+1, . . . , Om) = Yl "^i^i ~ ^j) II ^-(^i) II • • • , On+l] On, ■ ■ ■ , Ol) 

i<j i=l i=n+l 
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Due to integrability it factorizes into pairwise scatterings and individual transmissions: T_(^) 
and T+(— 6'). We parametrize T+ such a way that for its physical domain (9 < 0) its argument 
is always positive. Transmission factors satisfy unitarity and defect crossing symmetry p!3| 

T+i-e) = T:\e) T^{e) = T^{m - e) (3) 

The multiparticle transition amplitude can be derived by introducing the defect operator 
and the following relations in the ZF algebra: 

A+{e)D+ = T_{e)D+B+{9) ; D+B+^-O) = T+{e)A+{-e)D+ 

A defect is parity symmetric if T_(6') = T+(6'). Clearly A^{9 = 0) satisfies the properties of 
with T_(6') = S{9) = T+{9) . Thus a standing particle can be considered as the prototype 
of a parity symmetric defect. 

2.2 Coordinate dependence of the form factors 

The form factor of a local operator 0{x,t) is its matrix element between asymptotic states: 

where the adjoint state is defined to be 

C'^ • • • ■■■A\ = m{9\) ■ ..B{9\^^)DA{9\) . ..A{9[) 

Strictly speaking the form factor is defined only for initial/final states (i.e. for decreas- 
ingly/increasingly ordered arguments) but using the ZF algebra we can generalize them for 
any values and orders of the rapidities. 

The multiparticle asymptotic states are eigenstates of the conserved energy. This fact can 
be formulated in the language of the ZF algebra as 

[H,A+{9)] = mcosh9A+{9) ; [H,D+] = eDD+ 

In the second equation we supposed that the vacuum containing the defect has energy cd- 
Classical considerations together with the topological nature of the defect suggest the existence 
of a conserved momentum with properties 

[P,A+(e)] = msinh^A+(e) ; [P,D^]=pdD^ 

Thus, opposed to a general boundary theory, the defect breaks translation invariance by having 
a nonzero momentum eigenvalue Pd and not by destroying the existence of the momentum 
itself. As a consequence the time and space dependence of the form factor can be obtained as 

(6'^,, . . . , 6'^/+i; 9'^,, . . . , 9[\0{x,t)\9i, ...,9n] 6'„+i, ...,9m) = 

itAE-ixAP rpO ( nl fl' ■ f)' fi' \fi ft ■ ft ft \ 

^ ^ {n' ,m')(n,m)\'^ni+mn •••) nn •■■) C'lrl? • • • j ^n, Un+m) 

where /S.E = m{J2j cosh.9j — cosh.9j,) and AP = m(X]j sinh6'j — X)/ sinh 6'^/). The very 
same simple space and time dependence can be seen also in a theory without the defect and 
it is substantially different from what we would expect from a general boundary theory where 
even the one point function has a nontrivial space-dependence. These considerations remain 
valid for operators inserted at the defect 0(t) = O{0,t), too. 
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2.3 Crossing transformation of defect form factors 



The properties and analytical structure of the form factor F(n'm'),{nm) can be derived via the 
reduction formula from the correlation functions similarly to the boundary case [9]. Instead of 
going to the details of the calculation of [9] we note that all equations follow from the defining 
relations of the ZF algebra and the locality of the operator [0(0, 0), A+(6')] = except the 
crossing relation. This missing relation reads as 



and can be obtained as follows: We fold the system [13j to a boundary one: B^{6) ^ B^{—9), 
and consider A'^ and as creation operators of two different type of particles which scatter 
trivially on each other. Now we apply the crossing equation of for the resulting boundary 
form factor [9]. If we fold back the system to the original defect theory we obtain the defect 
crossing equation above. 

By analyzing the crossing equation of the particle ^4+ instead of B^ we obtain 



pO (n/ nl . a' 



(n/+l,m/){n— l,m) \^nl+mn 



On ; 9. 



., Or, 



Or,: 6. 



Or, 



This crossing equation can also be obtained from (l2|). Using any of the crossing equations above 
we can express all form factors in terms of the one-sided form factors: 



-^(n,m)(^l' 



-^(0,0)(n,m)l 



1^1, 



• ) ^ni Oj2-{-lj •••) On^rn 



on which we focus in the rest of the paper. The properties of this form factor follows from 
the ZF algebra relations and from the crossing relations and we postulate them in the next 
subsection as axioms. 

2.4 Defect form factor axioms 

The matrix elements of local operators satisfy the following axioms: 
I. Transmission: 



-^(n,m)(^l5 • 



On] On+1, On+m) — ^- (^n)-^(n-l,m+l) (^li ■ ■ ■ i On-l] On-, On+m) 





By means of this axiom we can express every form factor in terms of the elementary one 

Fni^^l^ ■■■,0n) = -^(n,0)(^l' • • • 7 ) 
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It satisfies the further axioms: 
II. Permutation: 

F^{9i, . . .9i, 9i+i, . . . ,9n) = S{9i — 9i+i)F^ {9i, . . . 6'i+i, 9i,..., 9n) 



III. Periodicity: 




The physical singularities can be formulated as follows. 
IV. Kinematical singularity: 




V. Dynamical bulk singularity: 

-iRes6)'=6iF^2(^' + — , ^ - — , ^1, . . . , = TF^^^{9, 9i,...,9n) 



-i Res 
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where T is the 3 particle on-shell coupUng. 
VI. Dynamical defect singularity: 

-iRese=iuF^+^{ei, ...,6^,6)= tgF^iOi, . . . , 



-i Res 




where g is the defect bound-state coupling. 

A few remarks are in order: Although the form factor axioms (H-V) are the same as the 
axioms of the form factors in a theory without the defect [I], the axioms (I, VI) are different 
and in general defect theories will have different solutions. An exception is the invisible defect 
T_{6) = 1 when we recover the usual form factor equation providing a consistency check for 
our axioms. Another consistency check can be obtained by considering a standing particle as 
the defect. Then T±{6) = S{9) and the two additional axioms become part of the old ones: 
(I, VI) will be special cases of (II, V), respectively. 

3 Form factor solutions, two point functions 

In this section we determine the solutions of the form factor equations for operators localized 
in the bulk and at the defect. For operators localized in the bulk the solutions can be built up 
form the bulk form factors and from the transmission factors. Using these form factor solutions 
we determine the spectral representation of the two point function for the situations when the 
operators are localized on the same or on the opposite sides of the defect. Finally, for operators 
localized on the defect we outline the strategy for the general solution. 

3.1 Bulk operators 

The form factor axioms for Fn are the same as in the bulk so we expect to use the bulk form 
factor solutions. Clearly we have to make a distinction whether the operator are localized on 
the left, or on the right of the defect. If the operator is localized on the left then particles 
arriving from the left can reach the operator without crossing the defect. Since the defect is 
topological we can change its location without altering the form factor (as far as we do not 
cross the insertion point of the operator). Shifting then the defect far away we expect to obtain 
the form factors of the bulk theory. 




7 



Thus we can conlude that for the initial state > . . . > > the defect form factor 
coincides with the bulk form factor. Let us denote the solutions of the bulk form factor equations 
by Bn{9i, . . . , On). Then we claim that for an operator localized on the left (C^ for short) we 
have 

Fni^l, ■■■,0n) = Bn(Oi, ...,6n) = F^fi){(^l, ■■■,&n',) (4) 

By using the transmission axiom and the crossing relation we can express all other matrix 
elements in terms of the bulk matrix element and the transmission factor. If the operator is 
localized on the right of the defect (O^) then, by similar argumentations, we expect the defect 
form factor to coincide with the bulk form factor for particles coming from the right. Those 
initial states have the ordering > 6*1 > . . . > 6'„ and the form factor is then 

-^(0,n)(; ^1) ■ ■ ■ ) ^n) = Bn{Oi, . . . ,0n) 

The solution for the elementary defect form factor for operators localized on the right thus 
turns out to be 

F>(^i, ...,en) = '[{T_{ei)Bn{eu ...,9^) (5) 

i 

which satisfies all the bulk form factor axioms but does not coincide with the bulk form factor 
solution. 

Having calculated the form factor solutions we use them to construct the two point fimctions, 
which, for operators localized on the opposite side of the defect, will be intrinsically different 
from the one without the defect. We analyze the following two point function 

(;|C»i(xi,ii)02(x2,^2)|;) 

where we denote by |;) the vacuum of the defect theory. Formally |;) = D~^\0). Now we 
insert the resolution of the identity. It can be composed both from initial and from final states 
and for definiteness we choose initial states. It is instructive to list the possible states. If 
we have no particles we have only the vacuum: |;). One particle states can be of two types, 
depending on whether the particle arrives from the left or from the right: |^;) for ^ > 
and 1;^) for ^ < 0. A general N — n + m particle state \9i, . . . ,9n;0n+i, ■ ■ ■ i^n+m) with 
^1 > . . . > 6'„ > > ^„+i > . . . > 9rn has to cover all possible cases ranging from n = to 
n = N . The two point function then can be written formally as 

oo 

(; \0,{x^,h)02{x2,h)\;)= |Oi(0, 0)|iV)(; |0*(o, 0)|iV)*e^^W(*-*^)-^(^)(— ^) 

N=0 

We have to specify the integration ranges for the multiparticle state N. Originally we have 
to integrate only for the multiparticle momentum range of the initial states. If we exchange 
the order from 9i > 92 to the nonphysical 92 > ^^ithen the form factor of Oi picks up a factor 
S{9i — 92) while that of O* the inverse factor S*{9i — 92), so the integrand is a symmetric 
function. For each integration with ^1 > we also have an analogous integration for ^1 < . 
Their contributions differ by a factor for the form factor of Oi and by the inverse T*{9i) 

for O2. As a consequence we can express the correlator in terms of the elementary form factors 
F„ as: 

{■,\Oi{xi,ti)02{x2,t2)\;) = 

d9i d9 



00 2 1"^ 



2tt I , 2-K 
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Although we transported the operators Oi and O2 into the origin, the form factor solutions 
remember whether the operators are localized on the left or on the right of the defect. 

If both operators are localized on the left, (xi < 0, X2 < 0) then the elementary form factors 
are the same as the bulk form factors ^ and we can conclude that the two point function is 
exactly the same as the bulk two point function 

(;|Oi(xi,ti)02(x2,t2)|;) = 

n=0 ^' J -00 J -oo 

This is intuitively clear: we can transport the defect to infinity without crossing any of the 
insertion points thus leaving invariant the two point function. When the defect is at infinity it 
does not infiuences the two point function which then has to be the same as in the bulk. The 
same result can be obtained when both operators are localized on the right of the defect. 

If the operators are localized on different sides of the defect (xi < 0, 0:2 > 0) then additinally 
to (lH) we also have to use ([5]). As a result the two point function is expressed in terms of the 
bulk form factors i?„ and the transmission matrix T_ {6) as 

{■,\Oi{xiM)02{x2M\;) = 

This is the main result of this section. This formula shows how the correlation function can 
be calculated in the presence of an integrable defect in terms of the transmission factor and 
the bulk form factors. It can be generalized to any correlators localized in the bulk using the 
resolution of the identity together with the exact form factors (111) and (JH]). It cannot be applied, 
however, for operators localized on the defect, which is the subject of the next subsection. 



3.2 Defect operators 

We have seen that although the minimal form factors are subject to the same requirement 
as the bulk form factors they are not necessarily the same. In this subsection we develop a 
general methodology to determine the defect form factors. Let us analyze them for increasing 
particle numbers: 

The first form factor is the vacuum expectation value of a defect field 

{■,m)\;)=Fo 
The one particle form factor is defined to be 

{■\0{t)\e;)=F^{e) ; {■\0{t)\;e)=T4d)F^{d) 

Contrary to the bulk case it has a nontrivial rapidity dependence: it is not natural to take Fi 
to be constant. In a parity invariant theory for a parity symmetric operators, for example, we 
have Fi{9) = T^{—6)Fi{—9). If parity is broken then Fi{9) can be an arbitrary defect condition 
dependent 27ri periodic function. The only restriction came from the defect bound-state axiom 
(VI): it must have a pole at iv whenever T_{9) has a pole correspondig to a bound-state. Let 
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us denote the minimal function which satisfies this requirement by d{6). The general form of 
the one particle form factor is then 



{■,\0{m) = die)Q?iu) ; u = e' 

where d{9) depends on the defect condition, while Qi{u) depends on the operator we are dealing 
with. 

The two particle form factor must also have a singularity at iu and additionally it satisfies 
the bulk form factor axioms so we expect it to be written into the form 



Ui + U2 

where fmin{S) is the minimal solution of the bulk two particle form factor equations 

fmin{0) = S{d)fmin{-d) ] fmin{iT^ ^ ^) = fmin{iTT + 0) 

Here Q2 is a symmetric function in Ui = containing a factor (M1+M2) to kill the denominator, 
what we introduced to conform with the kinematical singularity appearing at higher levels only. 
Taking into account the general paramterization of the bulk and boundary form factors together 
with the dynamical and kinematical singularity axioms we parametrize our minimal defect form 
factors as 

FniOu . . .,0n) = n ^^;~^^-^ g,(^,, ...,«„) (6) 

where Qn{ui, ■ ■ ■ ,Un) is a symmetric function expected to be a polynomial, if there is no bulk 
dynamical singularity. If there is such a singularity we have to include the corresponding 
singularity into fmin ■ The dependence on the defect condition is contained in d{9) with possible 
defect bound-state singularities, while the dependence on the operator in Q. If for instance the 
defect is the invisible defect with T± = 1 then d = 1 and we recover the solution of the bulk 
form factor equation as it should be. From the kinematical singularity equations recursion 
relations can be obtained among the polynomials Qn+2 and Qn- 

4 Model studies 

In this section we analyze the solutions of the defect form factor axioms for the free boson and 
for the Lee- Yang models. 

4.1 Free boson 

The purely transmitting free bosonic theory was analyzed in [16] as the limiting case of the 
sinh-Gordon theory. Let us recall its solution. The Lagrangian of the model reads as 



C = 0(-x) 
6{x) 



2^ ' ^ 2 " 



0(x) 



(cosh /i) + $^ ) + 2(sinh /i) 
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where $± are the fields living on the right/left part of the defect and is a free parameter. By 
varying the action we obtain the free equation of motion in the bulk 

{dl - dl)^^ = m^$± (7) 

and the defect condtions: 

±dt^±^dx^^ = m(sinh//)$± + m(cosh//)$=p (8) 

Since $± are free fields they have an expansion in terms of plane waves and creation/ annihilation 
operators 

where the a, a"*" operators are adjoint of each other with commutators: 

[a±{k), al{h!)\ = 27r2uj{k)5{k - k') 

They are not independent, the defect condition connects them as 

/ ,N ^ /,x / ,x ^ /,N msinh u ^ ia;(/c) 

a±{±k) = TUk)a^{±k) ; T^{k) = , V ; A; > 

m cosh ji — ik 

This shows that the defect is purely transmitting, that is we do not have any refiected wave. 
The transmission factor in the rapidity parametrization {k — rrici sinh 6) can be written also in 
the following form: 

T ie) = - ^^^Mf-f + f) ^ 1 + ^^ 

^ ' sinh(| + f + f ) \-wu 

where w = ie^ and u — e^. Sometimes we also use w = w^^ and u = u~^. In the next 
subsection we will set m = 1 and use dimensionless quantities. 



4.1.1 Form factors 

In the free boson model, we have the advantage that we can explicitly calculate the form factors 
of all the operators, and then check that they satisfy the defect form factor axioms. Additionally 
we can also confirm that we have as many polynomial solutions of the axioms as many local 
operators exist in the theory. 

We work with the Euclidien version of the theory {t = iy) and introduce complex coordinates 
2z — y + ix , 2z — y — ix . Using the explicit expressions of $±(-2, z) we are able to calculate 
the form factors. We analyze them for increasing particle number and compare to the form 
factor solution. 

The one particle form factors turn out to be: 

Ff- = (0|$_(2,^)|a_(^)) = e^"+^"" 

Ff+ = (0|*+(^,^)|a_(^)) = e""+"""r_(^) 

from which it is easy to calculate the defect form factors of the derivative of the elementary 
fields: 
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(0|9"$_(0)|a_(^)) = ; (0|9"$_(0)|a_(^)) = u 



(0|a"$+(0)|a_(^)) =m"T_(0) ; (O|a"<l>+(O)|a_(0)) = M"r_(0) 

where d = and d = 82. We can unify this notation by = 9". It is instructive to see 
how we can recover these form factors from the solution of the form factor axioms. Using the 
parametrization of the form factors in terms of d{9) see equation ((61), we know that at level 1 
the solutions of the form factor axioms have the form: 

F^i9) = die).Q,i9) 

Thus if we choose 



we obtain 

Qf *- (^) = _ yju) ; Qf *+ (6) = + wu) 

Naively it seems we have less polynomial solutions of the form factor equations than operators. 
We have extra relations among the form factors and they originate from 

dd^± = ^± ; -9<l>_ = + ; a$_ + = «}($+- $_) 

However, these relations are satisfied due to the bulk equation of motion and the defect 
conditions ([8]). Note that the form factor solutions are even more simple in terms of = 
and (f) = $+ — Actually is not idependent since 4> = 2wd(t). Their form factors read as: 

Qf <^ = ■ Qf ^ = nju-^^ 

Observe also that by changing the sign of the exponential coupling —e^ the left right 

fields are interchanged $± ^ $=p, as follows from the discrete symmetries of the Lagrangian. 
In the form factor boostrap the general parametrization without kinematical singularity is 

N 

Fn{9i, . . . , 6'Ar) = Jl d{9i)QN{xi, Xn) 
i=l 

Thus we can read off the corresponding form factor solution directly 

Qn = w'l^ • • • — wui) ... (1 — wuk){l + wuk+i) . . . (1 + wun) + permutations 

Since the scattering matrix in the free boson theory is trivial S = 1, the form factors of different 
levels are not connnected to each other, and in this way we solved the theory completely. 

In terms of the field the form factor solutions are exactly the same as in the bulk free 
bosonic theory: 

Qn = ■ ■ ■ u}f + permutations 

thus we obtain exactly the same number of polynomial solution of the form factor axioms as 
many independent local operator exists in the theory. 
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4.2 Defect scaling Lee- Yang model 



The scaling Lee- Yang model can be defined as a perturbation of the M.{2,5) conformal minimal 
model with central charge c = — ^. It contains two chiral representations of the Virasoro 
algebra, Vq, Vi with highest weights and — |, respectively. The fusion rules can be summarized 
as: iVQj = iV*Q = 1 and iV{^ = 1 for i = 0, 1 and all others are zero. The Hilbert space on the 
torus corresponds to the (diagonal) modular invariant partition function and contains modules 
corresponding to the Id and the ^{z,z) primary fields with weights (0,0) and (— |, — |): 

n = Vo(E)Vo + Vi®Vi (9) 



The only relevant perturbation by the field $ results in the simplest scattering theory with one 
neutral particle of mass m and scattering matrix [2T| 

^ ^ sinhg + isinf 
sinh 6* — i sin I 

The pole at ^ = y (with residue F^) shows that the particle can form a bound-state. The 
relation 

s{e + z^)S{e-z^) = s{e) 

however, implies that the bound-state is the original particle itself and the bulk bootstrap is 
closed. 



4.3 Integrable defects 

Two types of topological defects can be introduced in the M{2,5) minimal model [28l[27j. They 
can be considered as operators acting on the bulk Hilbert space ((91) commuting with the action 
of the left and right Virasoro generators. They have to act diagonally on each factor in ([9]) and 
satisfy a Cardy type condition. This leads to two choices which can be labelled by the same 
way as the bulk fields: (0,0) and (1, 1). After making a modular transformation the defect is 
inserted in space. The operators living on the defect can be described as [26] 

id 

For the topological defect labeled by (0, 0) the operator space turns out to be 

and coincides with the bulk Hilbert space. This defect is the trivial (invisible) defect. 
For the other defect labeled by (1, 1) we obtain 

= ® K) + V^i ® + K) ® ^1 + 2 VI ® Fi 

For each of the representation spaces we associate a primary field Jc?, f{z), ip{z), ^-{z,z), ^+{z,z) 
with highest weights (0, 0), (-1/5, 0), (0, -1/5), (-1/5, -1/5), (-1/5, -1/5), respectively. The 
nonchiral fields ^±{z,z) can be considered as the left/right limits of the bulk field ^{z,z) on 
the defect. 
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The bulk perturbation by $ in the defect conformal field theory does not break integrability. 
In the case of the trivial defect the transmission factor is simply the identity T = 1. In the case 
of the defect labeled by (1, 1) we can introduce a one parameter family of defect perturbations 
as well by properly harmonizing the coefficients of the (p{z), i^{z) and ^{z,z) terms. We plan 
to analyze this issue in a forthcoming publication. Related investigations with only defect 
perturbations can be found in [23]. The bulk and defect perturbed theory is integrable and can 
be solved by exploiting how the defect acts on integrable boundaries, see [16] for the details. 
To summarize, in the calculation the bootstrap relation 

T.ie + j)T.{e-j) = T4e) (lo) 

was used together with defect unitarity and defect crossing symmetry ((HI) to fix the transmission 
factor as 

sinhf- + ■>—) 

T4e) = [b + i][b-i] ; [x]=^ . ",^2 (11) 

smh(| + «ff - «f ) 

(Actually the inverse of the solution is also a solution but the two are related by the b ^ 6 + b 
transformation) . 

We also note that the defect with parameter 6 = 3 behaves as a standing particle both from 
the energy and from the scattering point of view. 



4.4 Defect form factors 

In this subsection we apply the general method developed in Section 3 to determine the form 
factors of the defect Lee- Yang model. The form factor can be written as 

F^{9„ e„) = n dm n (^;~^^-) g„(^„ ...,«„) (12) 

The minimal solution of the two particle form factor equation is well-known and reads as [2 

u + u^^ — 2 

fmin{0) = — ■ -— - v{i7!- - 0) v{-m + 6) 

u + u ^ + 1 



where 



v{e) =exp\2 —e 
[ Jo ^ 



dx wx sinh | sinh | sinh | 



sinh^ X 



We also included the pole corresponding to the dynamical singularity equation by the denomi- 
nator. We choose the normalization of the form factors as in the bulk 



H„ 



4^3 \2'2v{0) 



Qn{ui, • • • , Un) is expected to be a symmetric polynomial in Ui and Ui. 

Let us turn to the determination of d{9). Due to the defect dynamical singularity for 
Fn{9i, . . . ,9n) the defect dependent term d{9) must have a pole whenever T^{9) has a pole. 
Similar equation is valid for Fo^„(; 9i, . . . , 9n) = Y[Tl^{9i)Fn{9i, . . . , 6'„) at the defect bound- 
states poles of T^{9) . We will take into account that the transformation 6^6 — 6 exchanges 



14 



T^{6) with Tj^{6) and we expect that it acts in a similar way on the form factors. (For parity 
invariant operators). The minimal solution with these requirements turns out to be: 



^^^^ 4sinh(f + g(6-5))sinh(f + g(6-7)) + 2 cos(f - i^) ^/2> + uv + u-^v 

(13) 

where we introduced u — ^ and v — v ^. This function satisfies two relevant relations: 

d{e + m)d{e) 



1 - 2 cos(f - 2ie) 1 - v?y'^ - u~^u^ 
and 

Singularity axioms generate recursive relations between the polynomials. The kinemetical 
recursion relation is given by: 

Qn+2{-U, U, Ml, Un) = Dn{u, Ui, Mn)Qn(«l, •••,««) 

with 

u 



Y\^{UUJ + UiUj){uUJ — UiUj) — Y\^iuUJ — UiUj){uUJ + UiUj) 



where we introduced u — , iv — uj~^, while the bound state recursion relation is : 



n-l 



Qn+liuu;,UUJ,Ui,...,Un-l) = (uu + U Y[(u + Ui)Qn{u, Ui, Un-l) 



i=l 



Now let us try to solve these recursions. 
4.4.1 Solutions 

Since Qn{ui, ...,Un) is supposed to be a symmetric polynomial, it is useful to introduce the 
elementary symmetric polynomials a^\ui, ...,Un) which are defined through the generating 
function: 

n n 

n(H + ^^.) = E^""M"H^l'-'^n) 

i=l k=0 

By means of these functions the kinemetical recursive relation for Qn reads as: 

(-l)"+^(5„+2(-«, U, Ml, Mn) = {u'^l^^ - 1 + M"^Z/"^)M^i)„(M, Ml, ...M„)Qn(Ml, Mn) 

with 

k=l m=l,odd V 3 / 
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We are going to find the form factors of the operators ^±{z, z), ip{z), i^{z) and their descendants. 
We can choose ^± as the defect limits of the right/left bulk fields, thus we know already all of 
their form factors. Taking into account the explicit form of d{9) together with T_{9) we find 

Qi~ = vai + vai + ; ^ = vcri + v^i - 

For the two particle form factors we get 

Qt' = o-i{v'^a2 + V^vai + (Tiai + 1 + Vd,vai + 

and 

= cri(f^cr2 — V^fCTi + (TiCTi + 1 — V^vai + i''^(T2) 

where we used the solution of the bulk form-factor equation Qf = ai. They both satisfy 
the dynamical recursion relations. The asymptotics of the solutions for x ±00 refiect the 
dimensions of the fields (— —\)- 

In order to calculate the general form factor Qn of the fields we can take their bulk form 
factors 5* from [24l [25] and rewrite it into the form (fT2l l with the defect part given by ([Hj) . 
The resulting form factors Qn are all polynomials and have the right asymptotic properties. 
Similar considerations hold for the descendants of the identity operator. However, calculating 
the left and right limits they can differ. This is due to the defect perturbation and the relation 
between the two must be expressed in terms of the form factors of other fields which follows 
from the fact that in the UV these operators have the same limit. The relation between the 
left and right limits of the energy momentum tensor is the analogue of the defect condition we 
have seen already in the case of the free boson. 

Now we would like to describe the two chiral fields ip{z) and ip{z) which have dimensions 
(— i,0) and (0, — |). The corresponding solutions at level one with the right asymptotics have 
to form 

Qi =(ri ; Q{ = Ox 

They are related by the u ^ transformation. Using our recursion relations we find the 
related solutions at level 2 

= cri(fO-2 + z/) ; Q2 = ^l{l^^2 + v) 
Those solution are not unique as they might mix with the kernel solutions 

Ki = o-i(crio-i - 1), K2 = cTiicrl - 0-2), K3 = (7i{al - 02) 

and their descendants. Here we included the factor ai in all cases as F2 must not have a 
dynamical pole, so the denominator Xi + X2 has to be killed. Interestingly all these kernel 
solutions can be expressed in terms of ip,(f and the two descendants — $+), (9($_ — 

$+). This is a nontrivial statement and shows that the form factor solutions at level 2 having 
degree less then four are in one-to-one correspondence with the operator content of the theory. 

4.4.2 Parity Symmetry 

In this part we analyze how the parity transformation acts on the form factor solutions. The 
action of the parity operator P on the operators can be written as 
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Operator 


Qi 


Q2 




uai + z/cTi + -\/3 




$+ 

-r 




o"i(f^o"2 — x/SfCTi + (TiCTi + 1 — ^/Si'd'i + y^a^) 




0-1 


cri(w2 + z/) 




^1 


cri(i?a-2 + t;) 



Table 1: The form factor solutions of the primary fields up to level 2 



The corresponding action on the form factors is 

P{Q\OmO\ ) = (0|O^| ;-^) 
By explicit calculations we checked that$_and $+ are parity even with 

P(0|$_(0)|^; ) = (0|$+(0)| -,-6) 
while on the contrary, ip and ip are parity odd with 

p(o|^(o)|^; ) = -(01^(0)1 --e) 

These relations have to be confirmed in the Lagrangian framework. 

5 Boundary form factors via defects 

In this section we intend to illustrate how defects can be used to generate new boundary form 
factor solutions from old ones. The underlying fusing idea for the reflection matrices can be 
explained as follows: Suppose we place an integrable defect with transmission factor T_(6') 
in front of an integrable boundary with reflection factor R{9), which satisfies unitarity and 
boundary crossing unitarity: 

ITT ITT 

R{-e) = R-\e) ■ R{--e) = s{2e)R{- + e) 

If we fuse the defect to the boundary the composite boundary system will be integrable and 
will have refiection factor 

R{e) = T+{e)R{e)T4e) 

which, due to the defect unitarity and crossing equations, will satisfy boundary unitarity and 
crossing unitarity. This idea has been used to calculate the transmission factors from the 
already determined refiection factors R, R in the sinh-Gordon and Lee- Yang models in [16]. In 
contrast, here we would like to use the fusion idea to generate new form factor solutions from 
old ones. For this purpose we suppose that we determined already the boundary form factors 
F^{9i, . . . ,9n) of a boundary operator O. It satisfies, besides the singularity [9j axioms, the 
following requirements: 
permutation 

F^{9i, . . . ,9i, 6'i+i, . . . ,9n) = S{9i — 9i+i)F^{9i, . . . , 6'i+i, 9i, . . . , 9n) 
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reflection 

and crossing reflection 

F^{e,, 02,..., e^) = R{tn - (2m -e^,e2,..., 9n) 

We claim that the fused form factor 

n 

F^{9,, . . . , = n T.i9.)F^{9,, ...,9n) (14) 

i=l 

satisfies the boundary form factor axioms of the fused boundary corresponding to the reflection 
factor R. 

Let us analyze them one by one. Since the extra factor is symmetric in 9i the permutation 
axiom is trivially satisfied. To show the reflection property we use defect unitarity 

R{e) = T+{9)R{9)T^{9) = T^{-9)-^ R{9)T^{9) (15) 

while for the crossing reflection we use defect crossing symmetry: 

RiiT: -9)= T+{iTT - 9)R{t7T - 9)T^{m -9)= T^{9)R{iTT - 9)T^{2m - 9)'^ (16) 

Now multiplying both sides of the reflection and crossing reflection equation by Y{iT-{9i) and 
using (fTSl) and (fT6l) the claim follows. Similarly one can show the satisfaction of the singularity 
axioms [9]. 

By this method form factor solution of a given boundary can be used to generate form factor 
solutions for the fused boundary. It is practically useful if we can follow the indentiflcation of 
the operators under the fusion procedure. This is the case for example if the operator in the 
UV limit commutes with the defect. Say for example if in the Lee- Yang model we take the form 
factors of the operators of the identity module on the trivial boundary p9], then by the fusion 
procedure we can generate the form factors of the same module on the fused boundary, just 
by multiplying the original form factor solutions with the product of the transmission factors. 



6 Conclusion 

In the paper we initiated the form factor program for purely transmitting integrable defect 
theories. We restricted our interest for a single particle type, but the extension of the program 
for diagonal bulk scatterings and diagonal transmissions is straightforward (see [30| in the 
boundary case). We laid down axioms for the form factors of operators localized both in the 
bulk and also on the defect. We determined the solutions of the consitency requirements for 
bulk operators in terms of the bulk form factors together with the transmission matrix. These 
form factors determine the correlation functions of bulk operators, which we elaborated in 
details for the two point functions. In the case of defect operators we gave the general form of 
the solutions and explicitly calculated for the free boson and for some operator in the Lee- Yang 
model. We also described how the fusion method can be used to generate new form factor 
solutions from old ones. 

In the analysis of the Lee- Yang model we observed relations between the defect operators 
which should have the origin in defect conditions. The lack of the Lagrangian definition of the 
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model prevented us to analyze this question. In order to acheive this aim one has to analyze 
the simultaneous integrable defect and bulk perturbations of the defect Lee- Yang model using 
conformal perturbation theory and establish the relation between the bulk and defect couplings, 
which maintains integrability. This approach then can be used to derive defect conditions which 
will provide relations between fields living on the defect. In pushing forward this program one 
has to solve the defect Lee- Yang model first. The expilicit knowledge of the correlation functions 
together with the structure constants will make it possible to bridge the operators appearing 
in the form factor program to their UV counterparts by analyzing the short distance behaviour 
of the two point functions obtained from the explicit form factor solutions. 

Another interesting problem is to see that we have as many polynomial solution of the form 
factor equation as many local operators existing in the theory. We have seen this coincidence 
in the case of the free boson. In the case of the Lee- Yang model the nontrivial mixing between 
the left and right degrees of freedom and the various cancellation between the leading order 
scaling terms prevented us to perform this analysis. Possibly a more careful analysis along the 
line of [29| would clear up this point as well. 

We have analyzed the free boson and the Lee- Yang model sofar. The method, however, 
has a straightforward application for the sinh-Gordon model adopting ideas from the boundary 
form factor solutions [311 [32]. 

The defect form factors in the Lee- Yang model can be tested by extending them for finite 
volume and comparing to direct TCSA data. They also can be used to build up finite tem- 
perature defect correlation functions. These are direct generalizatons of the related boundary 
analysises developed in [331134]. 
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